The Real Solution to Scalar Field Equation in 5D Black String Space 
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After the nontrivial quantum parameters f2„ and quantum potentials Vn obtained in our previous 
research, the circumstance of a real scalar wave in the bulk is studied with the similar method of 
Brevik (2001). The equation of a massless scalar held is solved numerically under the boundary 
conditions near the inner horizon re and the outer horizon re Unlike the usual wave function in 
4D, quantum number n introduces a new functions ^^^in, whose potentials are higher and wider with 
bigger n. Using the tangent approximation, a full boundary value problem about the Schrodinger- 
like equation is solved. With a convenient replacement of the 5D continuous potential by square 
barrier, the reflection and transmission coefficients are obtained. If extra dimension does exist and 
is visible at the neighborhood of black holes, the unique wave function may say something to 

it. 
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I. INTRODUCTION 



_ In 1974 [jy] and 1975 Stephen Hawking published his analysis of the effects of gravitational collapse on quantum 

Q^, fields, and predicted that black holes are not perfect black, but radiate thermally and eventually explode, named 
^ Hawking radiation. Since then, many people have used various methods and techniques to research black hole through 
the particles radiating from it, such as the simple Klein-Gordon particles and Dirac particles (for some early works, 
see Damour and Ruffini [sl and Chandrasekhar Q respectively). Recently, as for the former scalar field, Higuchi et 
^ I al. and Grispino et al. p gave the scalar field solution outside a Schwarzschild black hole, Brady et al. [7|], Brevik 
^ . et al. [1] and Tian et al. Q studied the Schwarzschild-de Sitter case and Guo et al. [l3| made further studies in the 
' Reissner-Nordstrom-de Sitter case and so on. 
\^ I Regarding higher dimensional background, in order to avoid the interactions beyond any acceptable phenomeno- 
' logical limits, people assume that standard model fields (such as fermions, gauge bosons, Higgs fields) are confined 
on a (3 + 1) dimensional hypersurface (3-brane) without accessing along the transverse dimensions. The branes are 
• embedded in the higher dimensional space (bulk), in which only gravitons and scalar particles without charges could 
' propagate under standard model gauge group. In this paper, the massless scalar field is chosen. Meanwhile, there are 
: also many works (for a review with lar ge e xtra dimensions see Ref. TT'l) focusing on Hawking radiation such as Kanti, 
' Frolov and Harris et al. [13 [3 [3 |l5j [TGi] . in which they employed both analytical and numerical techniques to 
^ ■ calculate various particles' greybody factors and different energy emission rates on the brane and in the bulk outside a 
k>( ' (4 + n)-dimensional Schwarzschild black holes and (4 + ri)-dimensional Schwarzschild-de Sitter black holes. These are 
rS \ two kinds of small higher dimensional black holes with horizon rn ^ L, where L is the size of the extra dimensions. 
' Such small black holes can be treated safely as classical objects. 

That the world may have more than four dimensions is due to Kaluza (1921) [l3| and Klein (1926) [IS], who 
realized that a 5D manifold could be used to unify general relativity with Maxwell's theory of electromagnetism. 
After that, many people are more interested in the higher dimensional gravitational theory. Here, we consider the 
Space-Time-Matter (STM) theory presented by Wesson and co-workers [l9| [1^. This theory is distinguished from the 
classical Kaluza-Klein theory for a non-compact fifth dimension, the 4D source is reduced from an empty 5D manifold. 
Because of this, the STM theory is also called induced matter theory and the effective 4D matter is called induced 
matter. That is, in STM theory, the 5D manifold is Ricci-flat while the 4D hypersurface is curved by the 4D induced 
matter. Mathematically, this approach is supported by Campbell's theorem which states that any analytical solution 
of N-dimensional Einstein equations with a source can be locally embedded in an (N -I- l)-dimensional Ricci-flat 
manifold [21]. In the framework of STM, people studied many astrophysical implications such as Perihelion Problem 
[2^ , Kaluza-Klein Solitons [1^ [l^l , Black hole [1^ , Solar System Tests [2(|] and so on. So far as we know, except for 
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our previous work [27|, no one has researched the radiation of 5D black holes in STM theory before. 

This paper is organized as follows: In Section II, the 5D black string metric and the time-dependent radial equation 
about R^{r,t) are given. In section III, by a tortoise coordinate transformation, the radial equation becomes a 
Schrodinger-like equation. According to the boundary condition and the tangent approximation, a full numerical 
solution is presented. In section VI, with replacing the real potential barriers around black hole by square barriers, 
the reflection and transmission coefficients are naturally obtained. Section V is a conclusion. 

We adopt the signature (+, — , — , — , — ) and put h, c ,and G equal to unity. Lowercase Greek indices fi,i> . . . will be 
taken to run over 0, 1, 2, 3 as usual, while capital indices A, B, C, . . . run over all five coordinates (0, 1, 2, 3, 4). 



II. THE KLEIN-GORDON EQUATION IN THE 5D SCHWARZSCHILD-DE SITTER SOLUTION 

Within the framework of STM theory, a class of exact 5D solutions, presented by Mashhoon, Wesson and Liu [26 
[29I 1 1911 ■ describes a 5D black hole. The line element takes the form 



3 



fir)dt' 



1 



fir) 



{d9^ 



de 



In this metric 
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where ^ is the open non-compact extra dimension coordinate, A is the induced cosmological constant and M is 
the central mass. The part of this metric inside the square bracket is exactly the same line-element as the 4D 
Schwarzschild-de Sitter solution, which is bounded by two horizons — an inner horizon (black hole horizon) and an 
outer horizon (one may call this cosmological horizon). This metric ([1]) satisfies the 5D vacuum equation Rab = 
and, therefore, there is no cosmological constant when viewed from 5D. However, when viewed from 4D, there is 
an effective cosmological constant A. So one can actually treat this A as a parameter which comes from the fifth 
dimension. This solution has been studied in many works (30j (sij [3^ (33j focusing mainly on the induced constant 
A, the extra force and so on. 

We redefine the fifth dimension in this model. 



(3) 



Then we use (HI '--^ (l3|) to build up a RS type brane model in which one brane is at ?/ = 0, and the other brane is at 
y = Ui. Hence the fifth dimension becomes finite. It could be very small as RS I brane model (ssj or very large as 
RS II model [13] ■ It has been noted that the relevancybetween STM and brane- worlds theories and some embedding 
of 5D solutions to brane models are studied in [3I] [33| [111 [3§|. For the present brane model, when viewed from 
a (^ or y = constant) hypersurface, the 4D line-element represents exactly the Schwarzschild-de Sitter black hole. 
However, when viewed from 5D, the horizon does not form a 4D sphere — it looks like a black string lying along the 
fifth dimension. Usually, people call the solution to the 5D equation '^^^Gab = ^b^^^dAB (A5 is the 5D cosmological 
constant) as the 5D Schwarzschild-de Sitter solution. Therefore, to distinguish it, we call the solution ([1]) a black 
string, or more precisely, a 5D Ricci-flat Schwarzschild-de Sitter solution. 

In our previous paper [27| . we have introduced a massless scalar field to stabilizing this black string brane model. 
Considering a single mode of the scalar field, the wave function for this mode may reach the maximum value but 
keep smooth and finite at the brane. Hence, a steady standing wave is constructed. A suitable superposition of some 
of the quantized and continuous components of L[y), which is one component of scalar field, may provide a wave 
function which is very large at y = and drop rapidly for y 7^ 0. Naturally, a practical 3-brane is formed at the y = 
hypersurface. 

With the redefinition ([3]), the metric ((T|) can be rewritten as 
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where y is the new fifth dimension. Expression ([2]) can be recomposed as follows 

f{r) ^ ^{r - re){rc - r){r - ro). 



(4) 
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The singularity of the metric (|4]) is determined by /(r) = 0. Here we only consider the real solutions. The solutions 
to this equation are inner horizon r^, outer horizon and a negative solution Tq — —{r^ -\- Tc). The last one has no 
physical significance, and rc and re are given as 



r, 



2 

(6) 



^^ = _cosr?, 
e = ^ cos(120° - 77) 



where = \ arccos(— 3Ma/A) with 30° <rj< 60°. The real physical solutions are accepted only if A satisfy AAf^ < i 

Following the work [221, the massless scalar field <f> in the 5D black string space, satisfies the Klein-Gordon equation 

□$ = 0, (7) 

where □ = -^-^^ iVdS"^^ dx^) d'Alcmbertian operator. We suppose that the separable solutions to Eq. 

([7]) are of the form 

'i* = -y^-RUr, t)L{y)Yi^{e, (8) 

where R^{r,t) is the radial time-dependent function, Yi,-a{6,(f)) is the usual spherical harmonic function, and L{y) is 
the function about the fifth dimension. The differential equation of i?^ (t, r) is 

1 2 ( Rui\ d f 2'-'^^( RujW 9 . . Ru. 



-r 



/(r) dt"^ \ r 

where is a constant which is adopted to separate variables {t, r, 6, (f>). Meanwhile, in our previous study [2^ . we have 
solved successfully the evolution equation L{y) in branes model. Its norm \L{y)\'^, which represents the probability 
of finding the massless scalar particle, gets an extremum on the branes. According to the standing wave condition in 
the bulk, the spectrum of is broken into two parts: one is the continuous spectra below |A and the other is the 
discrete spectra above |A. The quantum parameter f2„ is 

3 

Vt ' 4 

where n = 1, 2, 3 • • • and j/i is the thickness of the bulk. 



^" = — ^ + 7A, (10) 



III. THE BOUNDARY VALUE PROBLEM IN THE BULK 



A. The Schrodinger-like Equation 



A more important fact is the aspect of radial direction. In Eq. ([9]) time variable can be eliminated by the Fourier 
component e^"^* via 

RUr,t)^^..in{r)e-'^\ (11) 
where the subscript n presents a new wave function unlike the usual 4D case "^^jI Eq. ([9]) can be rewritten as 



-f{r)^{f{r)^) + V{r) 
ar ar 



whose potential function is given by 

V{r) = f{r) 
Now we introduce the tortoise coordinate 



Idfjr) ^ l{l + l) 
r dr 
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The tortoise coordinate can be expressed with the surface gravity as follows 



1 

2M 



1, / r 
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where 



Explicitly, we have 



2 



ore 

So under the tortoise coordinate transformation (fT4|) . the radial equation p2|) can be written as 



-4A/2V(r) 



(15) 
(16) 

(17) 
(18) 
(19) 

(20) 



which is of the form of Schrodinger equation in quantum mechanics. Because there are two different coordinates 
— r and x in it, people usually call it Schrodinger- like equation. The incoming or outgoing particle flow between 
inner horizon Tg and outer horizon Tc is reflected and transmitted by the potential V{r). Substituting the quantum 
parameters r2„ (jlOp into Eq. (jl3p . the quantum potentials are obtained as follows 



K(r) = /(r) 



ldf{r) ^ + ^ n^TT^ 
r dr yf 



-A 



(21) 



This means that the potential V{r) ([2T|) determines the evolution of the field 'i'^iin- 



B. The numerical solution 

Near the horizons and Vc, the tortoise coordinate x — !■ ±oo. According to Eq. ([5]) and Eq. pT|) . the boundary 
conditions are as follows 



V{r,) = V{r,) = 0. 



So near the horizons, Eq. (1201) reduces to 



r ^2 



dx'^ 



(22) 



(23) 



where the potential vanishes. Obviously, its solutions are e^*^*^'^^. In this paper, taking into account only real field, 
we choose the solution [8|] 



'^ujin = cos(2Mwx) 



(24) 



as a boundary condition near the two horizons ([6]). 

Because there are two coordinates — radial coordinate r and tortoise coordinate x coupled in Eg. ipDl) . the source 
transformation expression ([TS]) is too complex to be solved. In order to solve Eq. (PO]) conveniently, it is necessary to 
use an approximate method for transition 



f{x) = f{x). 



(25) 



As far as we know, there are two methods used frequently — one is tangent approximation Q, and the other is 
polynomial approximation [ol]. The former method is convenient for theoretical analysis and is adopted here. The 
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FIG. 1: The first three potential barriers Vn{x) of 5D black string with n — 1 (solid), n — 2 (dashed) and n — 3 (dash-dot). We 
use I = 1, M = 1, A = 10""^ and yi = 10"^^^ (a very large 5th dimension). Meanwhile, the potential barrier of the corresponding 
4D solution {Q — 0) is plotted for comparison. 



latter one involves too many polynomials approximation and we do not consider it in this paper. To the parameter 
A, we can always find an appropriate approximate method from the above. For a widely separated horizons model, 
we adopt the same value A = 10"'^ as appeared in Q Q. Then putting A into the Eq. we find the inner horizon 
fe = 2.00268M and the outer horizon Tc = 53.7435M. We employ the useful tangent approximation 



r[x) = - arctan 





(26) 



in which b = 2.6/ {rc — r^,) and d — {vc + re)/2 — 3 are two parameters to linearly fit r and x. Thus, the potentials 
Vn{r) ([2T|) can be converted into Vn{x), which are plotted in Fig. [TJ Clearly those potentials are highly locafized near 
X = and fall off exponentially at both inner horizon (x) and outer horizon rc (x) . It is characteristic in our higher 
dimensional scenario that potentials are higher with bigger quantum number "n" . 

Because such approximation (p6)) does not allow very large we shorten the distance of x to [-10,180], and hence 
boundary condition (|24p can be rewritten as 

*t^;,i(-10) = cos(20Afw), *^,„(180) = cos(360Mtj). (27) 

By using Mathematica software, combining the Schrodinger-like equation (|20p and boundary conditions (j27p . we can 
solve Eq. (j20p numerically as a boundary value problem. The amplitude versus the tortoise coordinate is illustrated 
in Fig. [21 It can be seen that the solution '^ajin{x) is similar to a harmonic wave without considering the decay factor 
1/r in expression ([8]). Taking into account the real surrounding, we use the tortoise transformation (jl5[) and plot the 
amplitude versus r in Fig. [3] Since the penetrating power becomes weaker, the wave is sparser than classical case Q 
near the black horizon. From the view of penetration of a square barrier, the difference can be demonstrated more 
clearly in the next section. 



IV. THE REFLECTION AND TRANSMISSION 

We assume that the particle fiux with energy E bursts towards a square well along the positive direction of x axis, 
where the potential is 

F(x) = ( -i<^<-2, (28) 

^ ' 1 0, X < Xi OT X > X2- 



20 40 60 80 

X 

FIG. 2: Variation of the field amplitude versus tortoise coordinate x with M = 1, n=l, tj = 1, / = 1, A = 10~^ and j/i = 10^' 
(a very large 5th dimension). 




10 20 30 40 50 



FIG. 3; Variation of the field amplitude versus r with M = 1, n=l, cj = 1, Z = 1, A = 10 ^ and yi = IQ^^^ (a very large 5th 
dimension). The waves pile up near the outer horizon. 

From the view of quantum mechanics, considering the wave behavior of the particles, this process is similar to 
scattering on the surface of propagation medium with thickness of \x2 — xi\. Parts of them are transmitted and parts 
of them are reflected back. According to statistical interpretation of wave function, whether the energy E > V[) ox 
not, there is definite probabilities to transmit and reflect by the potential. The reflection and transmission coefficients 
denote the magnitude of those probabilities. 

As mentioned above, it is necessary to replace the continuously varying potential barrier with a discontinuous barrier 
of constant height in analytical work. Therefore, the usual refiection and transmission coefficients can be obtained. 
With the method of [l^l and We suppose a scalar wave propagates from —oo to +cx), which is illustrated in Figl6l 
The same denotation is cited here, namely associating "1" with the incoming wave in the region — oo < x < xi, "2" 
with the potential plateau xi < x < X2, and "3" with the outgoing wave in the region X2 < x < +oo. Hence, potential 
V{x) in Eq.dnni) reduces to 

{Vi, — oo < X < xi, 
V2, Xi<X< X2, (29) 
Va, X2 < X < +00. 

Following those square barriers, the solutions to the Eq. (|20|) are 

{^^gifciK _|_ 5^e-«fcia:^ _oo < X < Xi, 
a^e'k^^ + &2e-*'=^^ xi < x < xa, (30) 
^^gjfcax^ X2 <_ X < +00, 



where ki = •w4M^(a;^ — V-i) (i — 1,2,3) are the wave numbers, and bi are the undetermined coefficients to the 
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FIG. 4: The reflected log J? versus the height H (or V2)- The same width d = 14 with n = 2 is adopted. Here, M = 1, I = 1, 
A = 10~^, and yi = 10^''^. The resonant points are marked by the star symbol (*). The n=2 case is represented by a black 
point. 



10-^ 
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10 r 



2 4 6 8 10 12 14 16 18 20 

the width d 

FIG. 5: The reflected logi? versus the width d. The same width H = 0.1093 with n = 2 is adopted too. Here, M = 1, Z = 1, 
A = 10~^, and yi = 10^''^. The resonant points are marked by the star symbol (*). The n=2 case is represented by a black 
point. 



solutions . Then we define the reflection coeflacients for the plane interfaces dividing two media 

where Zij = are the real impedance ratios between medium i and j. So the width of the barrier is d = X2 — xi and 
the height of the square barrier is if = t^. So in this model the reflection coefficients R and transmission coefficients 
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-20 -10 10 20 30 40 50 

X 

FIG. 6: Replacement of the real 5D quantum potential barrier Vi (solid) by square barrier with M = 1, Z = 1, A = 10~'^, and 
j/i = 10'^''^ (a very large 5th dimension). The replacement of the corresponding 4D solution (S7 — 0) is also plotted by dashed 
line for comparison. 



T are given as ^ 



R 
T 



Rl2 + ^23 + 2y/Ri2R23 COs(2fc2d) 

1 + R12R23 + 2V-Rl2i?23 cos(2fc2d) ' 

1 16 



(1 + Zi2)2(l + Z23)2 1 + R12R23 + 2ViCT^COs(2fc2d) ' 



(32) 
(33) 



It is well known that in quantum mechanics |4]| . the resonant effect frequently occurs in the square barrier. After 
the particle flux enters into the square barrier, the particle flux is reflected and transmitted again both by the walls 
X2 and xi (Fig. [6]). If wavelength satisfies the resonant condition, the resonant tunnelling would be arisen here. 
As a matter of fact, the waves after scattering many times have the same phasic. Their coherent states are stacked 
together. Hence, the amplitude of transmission waves increase evidently. Naturally, the reflection is weak more. Also, 
one can read this feature directly from Eq. (j32p and Eq. (j33p . which contain cosine functions. So there is no surprise 
that Fig. [Hand Fig. \5\ which describe logi? vs. H and logi? vs. d respectively, take oscillating like forms. When 
logi? gets the extremum points (stars), the resonant tunnelling is arisen in the square barrier. At the same time, the 
case of n = 2 (two black points) is presented in the both figures, too. 

In order to get reflection and transmission coefhcients in this 5D black string, we have to find the approximative 
square potentials to replace the continuous potentials Vn- So the replacements is plotted in Fig. [S] For clearness of 
the figure, here we only consider the 4D (O = 0) and the first 5D quantum potential {n — 1). In the case of n = 1, 
we read V2 — 0.1000 as a reasonable value for the potential plateau in Fig. [6l and choose the incoming wave number 



TABLE L The reflection and transmission coefficients 



mode 




Xl 


X2 


d 


V2 


R 




T 


4D {n = 


0) 


-4 


6 


10 


0.0950 


2.4 X 10" 


'3 


0.9975 


5D (n = 


1) 


-4 


6 


10 


0.1000 


2.7 X 10" 


-3 


0.9973 


5D (n = 


2) 


-4 


10 


14 


0.1093 


2.5 X 10" 


-4 


0.9997 


5D (n = 


3) 


-4 


20 


24 


0.1267 


1.9 X 10" 


3 


0.9981 



1 See Ref [| for detail. 
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to be fci = 2 (Vi = 0) and fca = 2 (V3 = 0). For the case of 4D {il — 0), we employ the same width and read 
V2 — 0.0950. Meanwhile, using the same method, we can get the width and altitude of square barrier {n = 2) and 
(n = 3), respectively. In the end, every parameter in those modes is given in Table (|T|. Substituting those into Eqs. 
(|32| and ([33)) . the reflection and transmission coefficients {R, T) are obtained and listed in Table (jT| too. Because of 
the resonant effect, it is difficult to judge directly whether R (or T, notice R+T=l) would be smaller or larger with 
the change of n only by considering the width d or the height H. Careful contrast of the Fig. [4] and Fig [5] shows 
up some key differences that although they all oscillate integrality, the average order of R varies obviously with the 
height H rather than the width d. 



V. CONCLUSION 



In this paper we have solved the real scalar field ^'^/n and obtained reflection coefficients and transmission coefficients 
in the 5D black string space. We summarize what has been achieved, and make some further comments. 

1. The 5D black string solution presented by Mashhoon, Liu and Wesson in Ref. [1^ is exact in higher dimensional 
gravitional theory. In the work of paper [l^l, the Klein-Gordon equation has been successfully separated into three 
parts, <I> ~ R^{r^t)L{y)Yim{d , (f) , corresponding to the variables (r, f), the fifth coordinate y, and the usual spherical 
harmonic coordinates {0,ip). We notice that the extra dimension y or ^ affects the usual 4D by the parameter fl, 
which is contained in the potential p3p . The wave function L(jj) gets its extremum on two branes. Then by using 
the steady standing wave condition, the quantum potential Vn are obtained after quantizing the parameter f2. So the 
effective 5D potential (HI]), which is different from the usual 4D case, determines the evolution of massless particles 
around the black hole (note that when = 0, the potential (|13p reduce to the usual 4D Schwarzschild-de Sitter black 
hole potential 8]). Naturally, we can get a new radial wave function ^'^j/n being homologic with '^^i 0] in 4D. 

2. One purpose of this paper is to find the difference between '^^in and ^j^/. We study scalar particles scattering 
around our black string solution ([T]) in the bulk and get the scalar field solution under the tangent approximation. Eq. 

describes the evolution of one dimensional transmission of a wave through the potential barrier. According to 
Eq. effective potential (PTjl and boundary condition a full boundary value problem is presented. Because 

of the complex potential (|2ip and the approximative replacement (I26p , we only give the numerical solution instead of 
a analytical one. 

3. From the view of mathematics, an oscillatory cosine function cos(2fc2(i) or cos{2^yAM'^{oJ'^ — H)d) is contained 
both in R ([32]) and T ((33]) . It is obviously that the function R{H, d) and T(i?, d) are not monotonic function. We can 
see it clearly from Fig. [5] and Fig. [5] Considering the different replacement of square barrier obtained through the 
different quantum potentials, we only show the case n=2 (black points)in the figures. It is obvious that it is not easy 
to distinguish which one is smaller or bigger just by H or d. Since that our replacement is a approximate method 
and its width and height are adopted very roughly, it is not necessary to give this resonant condition. 

4. It is known that the hierarchy problem, which is why the characteristic scale of gravity [Mp ~ 10^ CeV) is 
about 16 orders of magnitude larger than the Electro- Weak scale (Wew ~ 1 TeV), is solved successfully by assuming 
the existence of extra dimensions both in ADD model [i^l [1^ [i3| and R-S model 34 1 [3^. It is considered that 
the extra dimensions would be really visible at the neighborhood of black holes and in the early epoch of universe . 
If it does so,, regarding the former, the extra dimensionality of space would emerge out. In the strong gravitational 
field — black hole or black string (this model), the best way to probe this is the Hawking radiation. It contains a 
lot of information such as the mass of black hole or black string, the topological structure of our brane-world, the 
magnitudes of effective cosmological constant A and so on. The appearance of wave function '^^in is purely a quantum 
effect of the fifth dimensions. In this paper we employ a approximative method to give the real scalar field solution 
in the bulk roughly. 

5. In this paper we only consider the quantum number n depended property of scalar field. Unlike the work p^ . 
we do not take into account the dependence on the effective cosmological constant A. We take A — constant instead. 
The reason is that in order to solve the Schodinger-like Eq. (PO]) . which has two coordinates variables r and x, we 
need to find a quasi function ([^5]) . Nevertheless, we have to choose the value of effective cosmological constant firstly 
in the key function (|15p . If the value of A varies, the approximate function has to be changed too. Furthermore, it is 
difficult to find out the approximate relation between the radial coordinate r and the tortoise coordinate x. So in this 
work, the previous value A = 10"'^ Q [13] and the tangent approximation (1^5]) is employed. Anyway, it is interesting 
to research this case and the further work is needed. 
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